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Abstract. We present UBVRI photometry of the open cluster NGC2422 (age ~ 10 s yr) down to a limiting 
magnitude V ~ 19. These data are used to derive the Luminosity and Mass Functions and to study the cluster 
spatial distribution. By considering the color-magnitude diagram data and adopting a representative cluster 
main sequence, we obtained a list of candidate cluster members based on a photometric criterion. Using a 
reference field region and an iterative procedure, a correction for contaminating field stars has been derived in 
order to obtain the Luminosity and the Mass Functions in the M = 0.4 — 3.5 Mq range. By fitting the spatial 
distribution, we infer that a non-negligible number of cluster stars lies outside our investigated region. We have 
estimated a correction to the Mass Function of the cluster in order to take into account the "missing" cluster 
stars. The Present Day Mass Function of NGC 2422 can be represented by a power-law of index a — 3.07 ± 0.08 
(rms) - the Salpeter Mass Function in this notation has index a — 2.35 - in the mass range 0.9 < M/Mq < 2.5. 
The index a and the total mass of the cluster are very similar to those of the Pleiades. 

Key words, open clusters and associations: individual: NGC 2422 - Stars: luminosity function, mass function 



1. Introduction 

' Accurate determinations of the stellar Initial Mass 
Function together with the star formation rate are funda- 
mental to understand star formation mechanisms and re- 
lated astrophysical problems. Since Salpeter's e stimation 
of th e IMF for stars in the solar-neighborhood l|Sah>eterl 
several investigations in Galactic and extra-galactic 
stellar systems seem to conver ge to a unive rsal IMF de- 
scribed by a broken power- law l|Scalolfl998|l . Many mod- 
els of stellar population, chemical evolution and galactic 
evolution adopt a priori a single IMF assuming its univer- 
sality, although thei r results are highl y sensitive to uncer- 
tainties in the IMF l|Kennicuttl . ll998|l . 

Testing the universality of IMF is a challenge for astro- 
physicists because several indirect evidences suggest that 
the IMF ought to systematically vary wit h the time du e 
to the different star forming conditions llLarsonl [l998). 
Nevertheless, no convincing proofs for a variable IMF still 
exist and evidences for uniformity of the IMF are deduced 
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by estimates f or different populations dScalol Il98fil Tl998: 
lKroupall200^l . 

However, a large scatter in the logarithmic power-law 
index, for stars more massive than 1 M@, is evident. In 
order to understand how large apparent IMF variations 
are due to uncertai nties inherent t o any observational es- 
timate of the IMF, iKroupa et all l)200l|l investigated the 
scatter, introduced by Poisson noise and dynamical evolu- 
tion of star clusters, of the power-law indices inferred for 
iV-body model populations. The resultant apparent varia- 
tion of the IMF defines a "fundamental limit" such that 
any true variation in the IMF that is smaller than this 
fundamental limit is not detectable. In addition, determi- 
nations of the power-law indices are subject to systematic 
errors arising from unresolved binaries. 

Being systems of coeval and equidistant stars with the 
same chemical composition, open clusters are key sam- 
ples in investigating the IMF and its possible spatial and 
temporal variations. Determination of the IMF of open 
clusters can however be challenging because of the con- 
tamination from background Galactic field stars. A fur- 
ther complication comes from the difficulty to transform 
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the observed Present Day Mass Function (PDMF) into the 
IMF using proper assumptions on the stellar and dynam- 
ical evolution, mainly affecting high and low-mass stars, 
respectively. 

In order to reduce these complications it is convenient 
to study young open clusters or star forming regions. In 
the low mass star range the IMF is more uncertain but 
stellar evolution effects are not important; in this mass 
range the Present Day Mass Function is representative of 
the IMF. 

Additionally, for the question of the universality or 
variability of the IMF, it is convenient to compare clusters 
with same age in order to highlight other possible para- 
metric dependence. For these reasons we choose to study 
the southern Galactic open cluster NGC 2422 which has 
an age comparable to the well studied Pleiades cluster. 

The equatorial (J2000.0) and galactic coordi nates of 
NGC 2422 are RA= 07 /l 36 m 36 s , Dec=-14°30' l|Lvngal 
Il987l) . /=230°97, 6=3°13, respectively. Estimates of the 
fundamental parameters of this cluster, such as age, dis- 
tance modulus and reddening we re given in the past b y 
various authors as summarized bv lBarbera et alJ l)2002f) . 

The age of ~ 10 s yr was est i mated using theoretical 
isochrones (jRoio Arellano et"aTl Il997|) : the most recent 
value of the distance d = 498+8g 5 P c ' corresponding to 
a distance modulus (M — m) = 8.48 + n '^, was deduced 
by H ipparcos measurements on 4 stars l|Robichon et all 
Il999|) while the most recent v alue of the r eddening 
E{B-V) = 0.088 was reported bv lDambisI l|l999h . 

Several papers have been devoted to NGC 2422 in the 
past : the most recent of them report Stromgren photome- 
try llShobbrookLll984HNisserl[l98llRoio Arellano et all 
Il997|) . while the available UBV photometric val- 
ues, either photographic and photoe l ectric , extend 
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only down to V ~ 16 (IZud. lumt 
Hoag et all Il96lt ISrnyth Nandvl 11962 
19661 llshmukhamedovl Il967|) . Mean photometric data 
and spectral classification from the former papers were 
compiled by iMermilliodl l|l986|) in a catalog of 212 ob- 
jects. Using this catal o g and 564 add i tional stars from 
llshmukhamedovl l)l967j) . lBarbera et alJ l)2002j) obtained a 
large literature-based compilation of measured data of 
stars in the field of NGC 2422; for some of these stars, 
X-ray counterparts were found. 

The layout of our paper is the following. We present, in 
Sect. 2, our photometry and astrometry and in Sect. 3 the 
method adopted to select the candidate cluster member 
sample. In Sect. 4 we describe how the Luminosity and 
Mass Functions of NGC 2422 and the spatial distribution 
of cluster members were obtained. Finally in Sect. 5 we 
summarize and discuss our results. 



2. Cluster UBVRI Photometry and Astrometry 

2.1. Observations and Data reduction 

The data used in this paper are CCD images in the UBVRI 
pass-bands collected at the 0.9-m telescope of the Cerro 



Tololo Inter-American Observatory (CTIO) on January 
30, 1997. The scale on the sky of the instrument is 2.028 
arcsec/pixel, for a total field of view of 1.15 x 1.15 square 
degrees, making such an instrument suitable to cover most 
of the apparent size of the cluster. The observing log is 
summarized in Table While the quality of seeing was 
limited, all images were collected in photometric condi- 
tions. Possible effects on the crowded field photometry, 
due to the limited seeing, have been ruled out by the ar- 
tificial star test described in Sect. 12.41 

Table 1. Log-book of the NGC 2422 CCD observations 



Filter 


Exp. Time 


Instrumental seeing 




w 


FWHM["] 


U 


30 


3.53 


B 


15 


3.24 


V 


15 


2.70 


R 


15 


2.63 


I 


23 


3.00 


I 


10 


2.99 



All the images were pre-processed in a standard way 
with IRAF, using the sets of bias and sky flat-field im- 
ages collected during the observing night. The instrumen- 
tal magnitudes and the positions of the stars for each 
frame were derived by profile-fitting phot ometry with th e 
package DAOPHOT II and ALLSTAR llStetsonl . Il987|) . 
Then we used ALLFRAME l|Stetsorlll994 to obtain the 
simultaneous PSF-fitting photometry of all the individ- 
ual frames. In order to obtain the transformation equa- 
tions relating the instrumental magnitudes to the stan- 
dard UBV (Johnson), i?/(Kron-Cousins) system, we also 
der ived the instrum ental profile- fitting photometry for the 
two lLandolll l|l992l) fields of standard stars SA95 and S A98 
observed during the same night. 

In order to obtain the total integrated instrumental 
magnitudes we derived aperture photometry for the same 
stars that we used to define the PSF, after the digital 
subtraction of neig hboring objec ts from the frames. We 
used DAOGROW ijStetsorl ll99C|) to obtain the aperture 
growth curves for each frame and to compute the aperture 
correction to the profile-fitting photometry. 

The transformation coefficients to the standard system 
were derived using transformation equations of the form: 

v = V + A + A 1 xX + A 2 x(V-I), 

b = B + B Q + Bi x X + B 2 x (B-V), 

i = I + C + C x x X + C 2 x ( V-I), (1) 

r =R + D + D 1 xX + D 2 x( V-R), 

u = U + E + E 1 xX + E 2 x( U-B). 

In these equations v, b, i, r and u are the aperture magni- 
tudes, already normalized to 1 sec exposure and X is the 
airmass. 

Due to the limited number of standard star obser- 
vations, we do not have a complete coverage of air- 
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Table 2. Coefficients of the calibration equations. 



2.3. The Color-Magnitude Diagram 



A 


Ai 


A 2 


4.637 ± 0.009 


0.14 


-0.012 ±0.010 


Bo 


Bi 


B 2 


5.574 ± 0.002 


0.25 


0.095 ± 0.005 


Co 


Ci 


C 2 


5.100 ±0.002 


0.05 


-0.008 ± 0.002 


D 


#1 


D 2 


4.352 ± 0.007 


0.10 


-0.029 ± 0.012 


E 


Ei 


E 2 


6.989 ± 0.002 


0.55 


-0.048 ± 0.008 



mass and we have not been able to derive the ex- 
tinction coefficients. In order to overcome this prob- 
lem, we used typical values for CTIO available at 
http:/ /www.noa o.edu/scope/ccdtime/ctio.db| both the 
adopted extinction coefficients and the best fit values for 
the zero points and the color coefficients are summarized 
in Table[3 Second order color terms were tried and turned 
out to be negligible in comparison to their uncertainties. 

The calibration process we have adopted is based on 
two sequential steps, following a procedure adopted by 
Stetson (private communication). First, 45 program stars 
on the NGC 2422 images were sel ected usi ng the condition 
(following a criterion described in lStetsorl 1993, Sect. 4.1) 
that each star has to be well separated from its neigh- 
bors, observed in all frames, and with a statistic index 
X, relative to the goodness of the PSF-fitting photome- 
try, less than 1.5. The photometric calibration based on 
the Landolt standards is applied to these selected stars 
only (which we refer to as local standards), which were 
then used (together with the Landolt standard stars) to 
calibrate the other program stars. 



2.2. Astrometry 

The astrometric solution has been computed using as ref- 
erence the recently released Guide Star catalog, Version 
2.2.01 (GSC 2.2). At beginning, our pixel coordinate list 
was matched to the celestial coordinate list of the GSC 
2.2 by projecting the celestial coordinates onto a plane 
and using as reference three stars for which we had both 
the pixel and the celestial coo rdinates from the Hipparcos 
catalog 1 l|Turon et allll993|) . 

An initial estimate for the linear transformation has 
been computed using the reference coordinates of the 1350 
matched stars. Then, a plate solution has been computed 
using the same matched pixel and celestial coordinates by 
fitting a power series polynomial (IRAF task CCMAP). The 
final accuracy is of 0.24 arcsec. Finally, the IRAF tasks 
CCSETWCS and SKYPIX were used to obtain the celestial 
coordinates of the total sample. 



In order to obtain the color-magnitude diagram (CMD) 
of the cluster, a selecti on based on th e sharp parameter 
was first done following lStetsonl l)l987[) . The sharp param- 
eter is related to the angular size of the astronomical ob- 
ject allowing to reject non-stellar objects as semi-resolved 
galaxies and unrecognized blended double stars, for which 
sharp is significantly greater than zero, or cosmic rays, for 
which sharp is significantly less than zero. We considered 
stellar objects those having the sharp parameter in the 
(—0.2 — ± 0.2) range. In this way we obtained a list of 
36 101 stellar objects for which we have V and / magni- 
tudes. For 35 732 of these we have VRI magnitudes, for 
35140 we have BVRI magnitudes and finally, for 33494 
we have UBVRI magnitudes. 

The V vs. V—I color-magnitude diagram for all the 
stellar objects in the NGC 2422 field is shown in Fig.[l](Ze/£ 
upper panel). Horizontal bars indicate the median errors 
in color, while vertical bars indicate the median errors in 
magnitude for bins of one magnitude. Clearly, the diagram 
is heavily contaminated by background and foreground 
stars, as expected from the position of the cluster in the 
Galactic disk. 

In order to determine the cluster main sequence mini- 
mizing the contamination effects, we show in Fig.{T\(right 
upper panel) the V vs. V—I color-magnitude diagram ap- 
proximately corresponding to the cluster core. The cen- 
troid of this circular region of ~ 16 arcmin of radius has 
been estimated using the stars with V < 11.5. 

We considered several theoretical isochrones avail- 



available at http://cdsweb.u-strasbg.fr/cgi-bin/VizieR 



able in literat ure llSiess et all EioTl iBaraffe et all Il99i 
iGirardi et all l2000i ICastellani et all Il999l) and we have 
found that the cluster main sequence is well fitted by 
the 100 Myr and solar met allicity (Z = 0.02 ) theoret- 
ical isochrone computed by ISiess et all <|000, SDF00), 
for highe r mass stars (M > .6Mm), and the one com- 
puted by iBaraffe et all l|l998l BCAH98) with a general 
mixing length parameter a = 1, for lower mass stars 
(M < 0.6-M©). Both isochrones were transformed into the 
V vs. V-I plane using the apparent distance modulus 
(m — M)v = 8.76, derived from the dista nce modulus 
(m - M) = 8.48 of iRobichon et all (|l999h and the in- 
terstellar absorption Ay = 0.273, from the stand ard re- 
lation R v = A V /E{B-V) = 3.1 (|Mathisl Il99rj(l where 
the reddening E( B-V) = 0.088 is the value calculated 
bv IPambisI |j999). The color excesses E( V-I) = 0.11, 
E(R-I) = 0.06, E(U-B) = 07 were derived us i ng th e 
relations for Ry =3.1 given bv lMunari fc Carrarol (|l996h . 
As shown in Fig. {left lower panel), the resulting theo- 
retical isochrone is in good agreement with the apparent 
cluster main sequence, confirming the literature cluster 
parameters, such as distance and reddening. 

Furthermore, using a set of theoretical isochrones by 
ISiess et all (|2000^ for different ages, and the well defined 
main sequence of bright stars in the B vs. U-B color- 
magnitude diagram, we verified that the NGC 2422 age 
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is ~ 10 yr, in agreem ent with the value given by 
iRoio Arellano et alJ tfl997^ . 

Finally, we matched our p h otome tric data with the X- 
ray sources of iBarbera et alJ l|2002F) . The matched stars 
are plotted with large squares in Fig. (right lower 
panel). There are two populations of X-ray sources: a blue 
faint population, dominated by background objects, and 
a group of sources belonging to the cluster main sequence. 
The photometric position of the X-ray members allow us 
to confirm our choice of the main sequence. 

In Fig. |2 we show the V vs. V-I, I vs. R-I, V vs. 
B-V and B vs. U-B color-magnitude diagrams for 
the stars in the NGC 2422 field. The curve in the V 
vs. V-I an d in the I vs. R -I color-magnitude dia- 
grams is the lSiess et alJ 1 2000|) theoretical isochron e ex- 
tended to lower stars using the iBaraffe et alJ (^98) the- 
oretical isochrone as described above. For the other two 
diagrams we only considered st ars with ctr < 0.0 3 and 
ajj < 0.03, respectively, and the ISiess et al.1 (|200C^ theo- 
retical isochrone because we use the V vs. B-V and the 
B vs. U-B color-magnitude diagrams to select only the 
bright cluster members. 



2.4. Data completeness and photometric errors 

Since we are interested in deriving the Luminosity 
Function from star counts and because of the limited see- 
ing of our data, particular attention has been devoted to 
estimate the accuracy of the photometry and the com- 
pleteness of the derived star list. For this, a list of artifi- 
cial stars was created and added to the original frames in 
order to compare the photometric results of the recovered 
artificial stars and the input values. 

In order to avoid overcro wding, the artificial star s were 
placed in a spatial grid (cf. iPiotto fc ZoccalA [l999) such 
that the separation of the centers in each star pair was 
four PSF radii plus 1 pixel. Using a random-number gen- 
erator for the V magnitudes, a list of 2470 stars (7% of 
the total sample) with a flat distribution of the instru- 
mental V magnitudes between 11 and 21.5 was created. 
These limits were chosen on the basis of the instrumental 
V vs. V-I color-magnitude diagram. The V magnitudes 
for each artificial star were converted to the other filters 
using the fiducial lines representing the V vs. V-I, V vs. 
B-V, I vs. R-I and B vs. U-B color-magnitude dia- 
grams. 

DAOPHOT's ADDSTAR routine was used to add these 
artificial stars into copies of the original data frames, with 
the appropriate frame-to- frame shifts in their position and 
brightness. Calibrated magnitudes were derived using the 
same photometric parameters and the same procedure de- 
scribed in Sect. 12. ll 

We estimated the completeness fraction as the ratio 
between the number of artificial stars recovered simulta- 
neously in the V, I and R filters and the number of added 
stars per one magnitude bin. This condition was imposed 
because the photometric selection of low-mass candidate 



cluster members was based on the use of V vs. V-I and 
I vs. R-I color-magnitude diagrams. 

We found that this ratio is equal to 1 down to V = 
18.26, while it decreases to 0.94 ±0.06 to the limit V ~ 19 
of our data. Therefore, we can conclude that our cata- 
log is complete down to V = 18.26 over the whole V vs. 
V-I and I vs. R-I color-magnitude diagrams while it is 
94% complete to V = 19, so that a ~ 6% correction will be 
required to the number of faintest stars when computing 
the Luminosity Function. 

Finally, photometric errors were estimated by the dif- 
ferences between the "observed" magnitudes and col- 
ors derived for the artificial stars and their known in- 
put values. We defined the external error <7 cxt for the 
V magnitudes and for the V-I and R-I colors as in 
IStetson fc Harris! l|l988j) ; the results are summarized in 
Table 

Table 3. External photometric errors from the artificial 
star photometry. 



AV 



Cext.V 0"ext,V-I <Text,R-X 



9.5 - 


10.5 


0.001 


0.001 


0.001 


10.5 - 


11.5 


0.001 


0.003 


0.003 


11.5 - 


12.5 


0.003 


0.004 


0.003 


12.5 - 


13.5 


0.006 


0.006 


0.004 


13.5 - 


14.5 


0.007 


0.009 


0.007 


14.5 - 


15.5 


0.013 


0.016 


0.010 


15.5 - 


16.5 


0.024 


0.025 


0.018 


16.5 - 


17.5 


0.043 


0.058 


0.039 


17.5 - 


18.5 


0.111 


0.105 


0.080 


18.5 - 


19.5 


0.225 


0.248 


0.098 



3. Photometric selection of Candidate Cluster 
Members 

In order to select possible NGC 2422 cluster members, we 
used only photometric informations as described in the 
following steps: 

1. In the V vs. V-I color-magnitude diagram, we se- 
lected as possible candidate members all those stars 
which, according to their uy and o~v-i errors, belong 
to a well defined strip in the CMD. The lower envelope 
of this strip follows the representative main sequence 
for the cluster (see Sect. l2.3|) while the upper envelope 
is displaced upward by 1 mag to include binaries. We 
have choosen 1 mag instead of the canonical 0.75 mag 
since our isochrone is, in some points of the diagram, 
slightly lower than the apparent "true" sequence of 
the cluster (see Fig. 1). With our larger strip we are 
sure to include X-ray detected stars that are probable 
members since significant X-ray emission is a common 
property of young stars. After this selection, our sam- 
ple contains 2059 possible candidate members. 
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1 2 3 1 2 

V-I V-I 



Fig. 1. The V vs. V-I color-magnitude diagram for the full catalog (left upper panel) and for stars within ^16 
arcmin from the centroid (right upper panel and lower panels). The curve is the adopted theoretical isochrone (see 
text) and the large squares indicate the positions of the stars which have X-ray counterparts in Barbera et al. (2002). 
Horizontal bars indicate the median errors in color, while vertical bars indicate the median errors in magnitude for 
bins of one magnitude. 



2. We considered at first only those stars (1895) for which 
UBVRI photometry was available. We constructed 
an analogous strip as above, in the / vs. R-I color- 
magnitude diagram and we rejected from our sample 
of initial possible candidate members 540 stars which, 
according to their 07 and o~r-i errors, do not belong 
to the strip in the / vs. R-I color-magnitude diagram. 

3. From the resulting sample, we considered only the 333 
stars having typical error in B magnitude less than 
0.03. We saw that this occurs for B < 16. Therefore, 



using the lSiess et"aH l|200f)() isochrone, we constructec . 
in the V vs. B-V color-magnitude diagram an analo- 
gous strip and we rejected those 58 stars that, accord- 
ing to their ay and <jb-v errors, do not belong to the 
strip. 

4. Subsequently, we allowed for the B vs. U-B color- 
magnitude diagram to select the bright stars with U < 
13.5 for which the errors in U magnitude are less than 
0.03 mag. Using the same criterion as above, we defined 
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Fig. 2. The V vs. V-I, I vs. R-I, V vs. B-V and B vs. U-B color-magnitude diagrams for the stars in the NGC 2422 
field. Horizontal bars indicate the median errors in color, while vertical bars indicate the median errors in magnitude 
for bins of one magnitude. The solid lines are the theoretical isochrones described in the text. In the V vs. B-V and 
B vs. U-B color-magnitude diagrams only stars with ub < 0.03 and o\j < 0.03, respectively, are indicated. 



a strip taking into account the isochrone's shape in this 
diagram thus further rejecting 195 stars. 
5. Finally, we considered the stars for which only VRI 
magnitudes are available. As above, we considered as 
possible cluster candidate members those belonging to 
the strip defined in the / vs. R-I color-magnitude 
diagram. 

In our final sample we retained only those stars with oi < 
0.3 mag. 

The selection procedure described above was chosen 
to take advantage of all the available photometric infor- 



mations and to avoid excluding too many possible cluster 
members. 

The photometric/astrometric catalog of the candidate 
cluster members, containing 1277 stars, is given in Table 
4 2 where we report RA and Dec (J2000.0) coordinates in 
decimal degrees, an identification number for each star, 
the U, B, V, R, I magnitudes and the associated uncer- 
tainties. 



2 available in the electronic format via the World Wide Web 
site http://cdsweb.u-strasbg.fr/ 
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Fig. shows these cluster candidates in the color- 
magnitude diagrams so far considered. Clearly, the sample 
includes contaminating objects that do not belong to the 
cluster. The approach to estimate and to deal with such a 
contamination will be described in the following section. 

4. Determination of Luminosity and Mass 
Functions and Spatial Distribution 

To construct the Luminosity and Mass Functions of the 
cluster we need to correct the luminosity distribution of 
our selected sample for field star contamination. To quan- 
tify the contamination we adopted an iterative procedure 
schematically summarized in the following steps: 

1. We defined as "field region" an area of our field of 
view where we expect to find only few cluster stars. 
We used such a region to quantify the field star con- 
tribution and then to obtain a first approximation of 
the cluster Luminosity and Mass Functions as well as 
of the cluster spatial distribution from stellar counts 
in a circular area centered on the c luster centro id. 

2. By fitting King's empirical profiles l|Kind .lT962) to the 
cluster spatial distribution we then estimated the total 
number of cluster members beyond our cluster region 
and the number of cluster members contaminating the 
"field region" . 

3. After correcting the field star luminosity distribution 
and the field star spatial distribution for cluster stars 
falling in the field region, we re-obtained the cluster 
Luminosity and Mass Functions and the cluster spatial 
distribution. By re-fitting King's empirical profiles to 
the new cluster spatial distribution we re-estimated 
the total number of cluster members beyond our field 
of view in order to correct the cluster Mass Function. 

4. We re-corrected the field star luminosity distribution 
and we repeated the above step 3. The entire proce- 
dure was rapidly convergent and indeed the "step 3" 
correction was applied only three times, since after the 
third application the previous and final corrected Mass 
Function differ by less than 1%, i.e. well within the in- 
trinsic uncertainties. 

4.1. First approximation of Luminosity and Mass 
Functions 

As summarized above, a first approximation to the con- 
tamination of foreground and background sources has 
been obtained using the lower region of our field where 
we expect to find only few cluster stars. In fact, as shown 
in Fig. 0] the cluster is mainly concentrated in the upper 
region of our field where we defined the " central cluster re- 
gion" corresponding to the circular area within the radius 
of 27.07 arcmin. A value for the cluster centroid position 
has been calculated as the median value of the X and 
Y coordinates of the bright stars (V < 12). This posi- 
tion corresponds to the equatorial coordinates (J2000.0) 
RA = 7 h 36 m 34 s .8, Dec = -14°29'18".5. We verified that 



the coordinates of the centroid do not change significantly 
if the V magnitude limit changes in the 11.5 < V < 13 
range. 




500 1000 1500 

X(pix) 



Fig. 4. A map of the brightest selected stars in the total 
field of view (1.15 x 1.15 square degrees). Biggest dots rep- 
resent stars with V < 9 while smallest dots represent stars 
with 15.3 < V < 15.7. The circle delimits the area chosen 
to estimate the Luminosity Function of the cluster, while 
the rectangle delimits the area chosen to estimate the con- 
tamination rate of foreground and background stars. 

We consider the rectangular region in Fig.0J of 67.6 x 
13.52 arcmin 2 size (2000 x 400 pixels), as representative 
of the Galactic field (the "field region"), i.e. we suppose 
that the observed luminosity distribution in the rectangu- 
lar region is given by 

MV) = N {fi (V) + N cl (V) ~ Nt, {V), (2) 

where Nf,o(V) is the "true" field star luminosity distribu- 
tion, while iV c i(V) is the cluster star contribution within 
the rectangular region. Therefore, the field star luminosity 
distribution in the "central cluster region" is given by 

N i (V) = N i (V)j, (3) 

where A and A are the area of the " central cluster region" 
and of the "field region" , respectively. 

The NGC 2422 Luminosity Function, ^(V), has 
hence been obtained subtracting the field star luminos- 
ity distribution, iVf(V), from the luminosity distribution 
of all the stars in our sample within the "central cluster 
region", N t (V). Fig. shows the contaminated luminos- 
ity distribution of the candidate members (dotted line), 
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Fig. 3. The color-magnitude diagrams for the possible photometric candidate members of NGC 2422. The sample 
has been obtained as described in Sect. EH 



the luminosity distribution of the field stars (dashed line) 
and the Luminosity Function of NGC 2422 obtained as de- 
scribed above and corrected for incompleteness as derived 
in Sect. 12.41 fsolid line). We found that the contamination 
rate from field stars is negligible for the brightest stars, 
but starts to become significant at V ~ 13 where it is of 
the order of the 40%. At fainter magnitudes this value in- 
creases with a peak of the order of ~ 90% around V ~ 16 
where the contribution of the field giant stars is domi- 
nant. The total number of the cluster members estimated 
at this step is 347. The Luminosity Function obtained at 
this step has been transformed into a first approximation 
of the Mass Function, £°(M), using the mass- visual mag- 
nitude relation derived from the same models that we used 



to fit the cluster color-magnitude diagrams. The resultant 
first approximation of the Present Day Mass Function, 
extending from 0.5 to 3.5 M Q , is shown in Fig.EJ 



As already mentioned, due to dynamical evolution, the 
angular size of the cluster is almost certainly larger than 
the investigated "central cluster region" hence this Mass 
Function is not representative of the whole cluster mass 
distribution. In order to estimate the correction to the 
Mass Function due to the dynamical evolution, we studied 
the spatial distribution of NGC 2422 as described in the 
following Sections. 
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Fig. 5. Luminosity distribution of the candidate mem- 
bers without correction for contamination (dotted line), 
the luminosity distribution of the field stars (dashed line) 
and the Luminosity Function of NGC 2422 (solid line) ob- 
tained as difference between the last two and corrected for 
incompleteness, as described in Sect. 12.41 
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Fig. 6. First approximation of the NGC 2422 Present Day 
Mass Function obtained in our central cluster region and 
without any correction. £°(M) values are given in number 
per logarithmic mass unit. 



4.2. First approximation of the Spatial Distribution 

The spatial distribution of the cluster members was stud- 
ied in order to investigate the cluster dynamical evo- 
lution and to determine the cluster dynamical parame- 
ters. As first approximation, we corrected the radial sur- 
face density distribution of candidate cluster members for 
foreground and background contaminating stars assuming 
that the field stars are uniformly distributed, i.e. 



dN „ 
— = 2pnr, 
dr 



(4) 



where r is the distance from the cluster center and p is the 
star distribution for unit of area in the rectangular region 
considered in Sect. 14.11 

In order to take into account the dynamical evolution 
and mas s segregation e ffects due to the energy equiparti- 
tion f see iKroupa et all 2001 and relative references), we 
subdivided the stars into 4 bins of V magnitude as in 
Table In Fig. we present the cluster surface density 
profiles in stars/arcmin 2 as a function of radius in parsec, 
assum ing the distance value of 497.5 pc bv lRobic hon et all 
l|l999l) . We fitted King's empirical profiles ijKinel Il962h . 



a 0.3 



8.5<V<1 1.125 




13.75<V<16 375 




0.3 



11.125<V<13.75 




a 0.3 



0.2 



16 375<V<19 




Fig. 7. The radial distribution of stellar density in 
NGC 2422, for 4 different bins of magnitude. 



given by 
p(r) 



VI + (r/rv) 2 \A + (rtAc) 2 



(5) 



to the surface density profiles. In this equation, p$ is a 
normalization constant, r c is the "core radius" determined 
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by the internal energy of the system and r t is the tidal 
radius where the cluster disappears. 

Th e tidal radius was calculated as in Ijeffries et alJ 
(|200l[) as r t = 1.46 Af c 1/3 , where M c is the cluster mass 
in solar masses. In order to estimate the total mass of the 
cluster we integrated the Mass Function obtained in Sect. 
14. H and found a value of M c — 357±35 M@, corresponding 
to the tidal radius value r t = 10.36 ± 0.34 pc. 

The results of King's empirical profile fitting, with the 
l-o" uncertainty estimates for the parameters, are reported 
in Tableland are shown in Fig. [7J Comparing the x 2 val- 
ues with the number of degrees of freedom (y — 6), we 
conclude that the fits are all acceptable. From the result- 
ing distribution we can see that, while the distribution of 
the brightest stars vanishes to zero as the radius increases 
(at the top on the left of Fig. [7J, the other profiles tend 
to non negligible values when fainter V magnitudes are 
considered. This result suggests an evidence of mass seg- 
regation of high mass cluster stars toward the center of the 
cluster and low mass cluster stars out of our investigation 
region of 27 arcmin of radius. This conclusion is confirmed 
by the increase of the core radius as the mass decreases. 

In order to estimate the fraction of low-mass stars be- 
yond the " central cluster region" , we extrapolated the sur- 
face density profiles as far as the tidal radius . Using the 
integral of equation El as given bv iKind l|l962l) . we calcu- 
lated, for each bin of magnitude, the total number of clus- 
ter stars within the tidal radius, indicated in Table [S] by 
TVtid- We compared the calculated total number 7V t id with 
the number of cluster members iV s c ^ r , detected in our sur- 
vey, that is within r sur = 3.7 pc and we found that while 
the stars with V < 11 are all in the investigated region, 
the fraction of the cluster stars lying outside our survey 
is 8% for 11.13 < V < 13.75, 22% for 13.75 < V < 16.38 
and 38% for 16.38 < V < 19.00. In order to verify the 
consistency of these numbers we also integrated King's 
empirical profile within the radius r sur , indicated by iVj™, 
finding consistent values with the number of the cluster 
stars found in the central cluster region (N^). 

The results of the King's empirical profile integration 
also suggest that the cluster star contribution within the 
adopted "field region", N C \(V) is non-negligible and that 
the estimated field star luminosity and spatial distribu- 
tions have to be corrected for the cluster star contamina- 
tion. In order to estimate the correction factor, we split 
the " field region" in a grid of 5 x 20 subregions and we cal- 
culated the King's profile surface density, p(ri), where 
is the radius of the centers of each subregion i. Therefore, 
the total number of cluster stars within the "field region" , 
per magnitude bin, is given by 

N cl (V)=A s J2p(ri), (6) 

i 

where A s is the area of each subregion. 

Finally, using equation J2J , we estimated a more accu- 
rate field star distribution as 

N it0 (V)=N i (V)-A s ^ i p(r i ). (7) 



4.3. Luminosity and Mass Functions 

Using the value of Nf,o(V) from Eq. JJJ) and the itera- 
tive approach described in Sect. we recomputed the 
Luminosity Function of the cluster taking into account 
the presence of cluster stars in the "field region" . 

This final Luminosity Function has been converted into 
a Mass Function that has been integrated to compute a 
new value for the total mass (M c ~ 385 Af ), correspond- 
ing to r t ~ 10.61 pc. Dynamical parameters were also 
obtained in the various iterations, and by integration of 
King's profiles we interpolated the ratio R = N^/Ntid, 
deriving a factor 1/R for the Mass Function to correct for 
the cluster stars lying beyond the investigated cluster re- 
gion. The corrected Mass Function has been computed by 
the expression: 

e c (M)=e(M)/R. (8) 

In Fig. [S]we compare the final Mass Function, £ C (M), 
with the Mass Function derived from the central cluster 
region but without correction for dynamical evolution. We 
also compare t hem with the Pleiad es Mass Function de- 
rived using the iLee fc Sunel lll995h Luminosity Function 
and the same mass-luminosity relation used for NGC 2422 
in the present work. 

2 1 0.6 0.2 0.1 M/M s 



CD 




a= 3.09 10. OH for 0.9<M/M e <2.5 
A NGC 2422 final 
O Pleiades 



0.5 -0.5 -1 

log (M/M Q ) 

Fig. 8. Comparison of the NGC 2422 Mass Function cor- 
rected for dynamical evolution effects (filled points) with 
the NGC 2422 Mass Function without correction for dy- 
namical evolution (empty triangles) and with the Pleiades 
one (empty circles). The line is the power law fitting the 
MF in the 0.9 — 2.5 M Q range. £(M) values are given in 
number per logarithmic mass unit. 
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Table 5. Results of King's empirical profile fitting to the radial surface density distribution. Columns 1 and 2 give the 
V magnitude ranges and the corresponding mass ranges; Col. 3 the fitted normalization constant, Col. 4 the assumed 
tidal radius, Col. 5 the fitted core radius, Col. 6 the % 2 of the fitting. Finally, Col. 7 lists the number of counted cluster 
members in our investigated region for each magnitude subset, while Col. 8 and 9 are the total number of cluster 
members predicted by the fitted model within the investigated radius (~ 3.7 parsec) and the tidal radius, respectively. 
The values are given for each of the four step of the adopted iterative procedure. 



AV 


M 


Po 


n 


r c 


x 1 


N c[ 

J v sur 


J v sur 


iVtid 






(M ) 


(stars/arcmin 2 ) 


(pc) 


(PC) 










ITERATION 


8.50 - 


11.13 


3.7 - 


1.7 


0.21 ±0.10 


10.36 


0.59 ± 0.22 


4.2 


34 


29 


33 


11.13 - 


13.75 


1.7 - 


1.0 


0.22 ±0.13 


10.36 


1.34 ±0.46 


8.5 


88 


76 


95 


13.75 - 


16.38 


1.0 - 


0.7 


0.27 ±0.24 


10.36 


1.49 ±0.73 


1.3 


106 


101 


129 


16.38 - 


19.00 


0.7 - 


0.4 


0.22 ±0.26 


10.36 


2.66 ± 1.48 


2.5 


119 


111 


164 


ITERATION 1 


8.50 - 


11.13 


3.7 - 


1.7 


0.21 ±0.10 


10.57 


0.62 ± 0.23 


4.0 


35 


30 


34 


11.13 - 


13.75 


1.7 - 


1.0 


0.23 ±0.13 


10.57 


1.36 ±0.44 


8.7 


93 


81 


103 


13.75 - 


16.38 


1.0 - 


0.7 


0.27 ±0.24 


10.57 


1.55 ±0.74 


1.3 


113 


108 


140 


16.38 - 


19.00 


0.7 - 


0.4 


0.24 ±0.25 


10.57 


2.65 ± 1.30 


2.5 


134 


124 


185 


ITERATION 2 


8.50 - 


11.13 


3.7 - 


1.7 


0.21 ±0.10 


10.60 


0.62 ± 0.23 


3.9 


35 


30 


35 


11.13 - 


13.75 


1.7 - 


1.0 


0.22 ±0.12 


10.60 


1.40 ±0.45 


8.7 


94 


80 


102 


13.75 - 


16.38 


1.0 - 


0.7 


0.27 ±0.24 


10.60 


1.58 ±0.76 


1.3 


115 


109 


142 


16.38 - 


19.00 


0.7 - 


0.4 


0.24 ±0.26 


10.60 


2.78 ± 1.41 


2.4 


136 


125 


191 


ITERATION 3 


8.50 - 


11.13 


3.7 - 


1.7 


0.21 ±0.10 


10.61 


0.62 ± 0.23 


3.9 


35 


30 


35 


11.13 - 


13.75 


1.7 - 


1.0 


0.22 ±0.12 


10.61 


1.39 ±0.47 


8.7 


94 


80 


102 


13.75 - 


16.38 


1.0 - 


0.7 


0.27 ±0.24 


10.61 


1.58 ±0.75 


1.3 


115 


109 


142 


16.38 - 


19.00 


0.7 - 


0.4 


0.24 ± 0.27 


10.61 


2.86 ± 1.47 


2.4 


137 


127 


195 



Ac cording to the recently proposed analytical IMF 
forms l)Krour>al [2002). we considered a multi-part power- 
law IMF to fit our data. We found that in the mass range 
0.9 < M/Mq < 2.5 the value of a is 3.09 ± 0.09 (rms) 
and 2.63 ± 0.16 for NGC 2422 and Pleiades, respectively. 
T he latest value is very close to th e value 2.67 derived 
bv lBarrado v Navascues et all l)200l[) for the Pleiades, us- 
ing ^tlelHa^I^^rEiil^dat a in the 0.6 - 10 M 
mass range. The two values here derived are in better 



agreement, within the errors, with thelScalol ( 199S ) value 

erj I 



a = 2.7 ± 0.3 rather than with the ISalpeterl ( 195E ) value 
a = 2.35. 

Using the final Mass Function, we also estimated the 
corrected cluster total mass (M C)Cor ~ 453 Mq) and we 
compared this value with the NGC 2422 total mass esti- 
mated within the "central cluster region", M c ~ 385 Mq. 
In the same mass range (0.5 — 3.5 Mq), the total mass 
of the Pleiades is Mpi = 467 ± 26 Mq , a value very simi- 
lar to the one of NGC 2422. We note that our estimate of 
the total mass of NGC 2422 is, in any case, a lower limit, 
since we do not have taken into account the presence of 
companions in photometric binaries. 



5. Summary and Conclusions 

Using UBVRI images, covering a field of view of 1.15x1. 15 
square degrees, we extracted a photometric and astromet- 
ric catalog of the stars in the field of NGC 2422 down to 



V ~ 19 (M ~ 0.4 Mq). Adopting a representative cluster 
main sequence to the color-magnitude diagrams, a photo- 
metric criterion was defined to obtain a list of candidate 
cluster members. A test with artificial stars allowed us to 
verify that our data are complete down to V = 18.26 while 
a correction of ~ 6% was necessary for the lowest V range 
(18.26 < V < 19.26) in deriving the Luminosity Function. 

We have defined a " central cluster region" within a ra- 
dius of 27.07 arcmin from the cluster ccntroid and a "field 
region" of 1.12 x 0.23 square degrees to estimate the con- 
tamination of background and foreground stars. An initial 
determination of the Luminosity and Mass Functions has 
been obtained assuming that all the cluster stars lie in the 
" central cluster region" . By applying an iterative proce- 
dure we have estimated the number of cluster stars within 
the "field region" in order to obtain a more accurate cor- 
rection of the cluster Luminosity Function for field star 
contamination. 

Evidence for mass segregation and energy equiparti- 
tion have been found from the spat ial distribution of the 
stars. Extrapolation of the Kins ( 1962) empirical model al- 
lows us to infer that, while all the stars with M > 1.7 Mq 
are inside the "central cluster region", a non negligible 
fraction of lower mass stars lies outside. In particular, by 
integrating the spatial distribution within the tidal radius 
we are able to estimate the number of cluster stars lying 
beyond the "central cluster region" and within the tidal 
radius. Therefore, we are able to estimate a correction to 
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the Mass Function in order to take into account the cluster 
dynamical evolution. 

The corrected Present Day Mass Function was com- 
pared with that of the Pleiades which is known down 
to the brown dwarf limit (0.08 M Q ). We found that, in 
the mass range 0.9 < M/Mq < 2.5, the Mass Function 
of NGC 2422 can be represented by a power law of in- 
dex a — 3.09 ± 0.08 (rms) comparable with the index 
a = 2.63 ± 0.15 obtained for the Pleiades in the same 
mass range. The index is also con sistent with th e data 
presented in the logM vs. a plot of lKroupal l|2002() . 

By taking into account the correction due to the dy- 
namical evolution, we computed a lower limit of the total 
mass of NGC 2422 as M c , cor ~ 453M Q . This value is sim- 
ilar to that of the Pleiades. 

A future spectroscopic study and a deep survey on 
a wider field of this cluster will allow us to find an in- 
dependent membership criterion and to extend the Mass 
Function below the mass limit of the present survey. 
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